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ABSTRACT 


The nutational stability of a dual-spin satellite under the influence of an applied 
reaction torque is analyzed. It is assumed that since the solar paddle attach- 
ments to the hub of the spacecraft are not rigidly locked, the effect of the solar 
panels can be replaced by a constant reaction torque acting on the hub of the 
spacecraft. This could result in the satellite having an equilibrium motion 
about an axis displaced from the nominal axis of symmetry. The variational 
equations of motion are developed about such an equilibrium position using the 
SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the 
main body is included. This nonautonomous set of differential equations are 
linearized and transformed to an autonomous set using the Lyapunov Reducibility 
Theorem. The stability of the kinematically similar system is examined numeri- 
cally using representative SAS-A parameters for the case when either pair of 
solar panels is assumed to be loosely attached. Stability is verified for small 
system nutation angles (0.1 degree) although the time constant associated with 
the least damped mode is approximately one order of magnitude larger than that 
for the SAS-A system without the influence of the reaction torques. 
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PRECEDING PAGE BLANK NOT FILMED NOMENCLATURE 


A, B, C 
A, B, C 
A (t) 


b. 

1 

b.' 

1 

Pl 

D 


H 

Iu 




di 


K 


Kr 

k 


k R 


L 


di 


= main body moments of inertia about the x, y, z axes respect 
= composite moments of inertia about the x, y, z axes respect 

= system matrix, appears when linear equations are written i 
vector form 

= unit vectors along the x, y, z axes respectively (i - 1, 2, 3) 

- unit vectors fixed to the nominal plane of the undeflected wb 
rotating with it 

= coefficients occurring in the linear equations 
= system matrix appearing in the kinematically similar syste 
equations 

= magnitude of system angular momentum vector, H 
= moment of inertia of main body about the bj axis 
= moment of inertia of rotor about its b. axis 
= moment of inertia of the pendulous damper about the b . axi 
= the restoring spring constant of the damper torsion wire su 
= the restoring spring constant of the rotor 
= the nutation damper damping (rate) constant 
= rotor damping constant 

= the applied reaction torque about the U axis 
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the damping torques about the b. axis 
height of damper plane above x, z plane 
the mass of the main satellite and the rotor 
the total system mass 
the pendulum end mass 
the mass of the rotor 

radial coordinate of a differential mass on the rotor 

the distance from the nominal spin (y) axis to the pendulum hinge 
point 

the length of the pendulum 

bounded nonsingular transformation matrix used to transform 
nonautonomous system of equations to an autonomous set 

spin rate of rotor relative to main body 

time 

principal axes of main satellite 

displacement coordinate of dm R from plane of undeflected rotor 


disc 



a = angle between the main body symmetry axis (b 2 ) and the main body 

angular velocity vector 

a a x = rotor deflection angles about the b 3 ',b 1 ’, axes respectively 

/3 v fi 2 = transformed rotor deflection coordinates in the kinematically 
similar system 


8, e, i// = variational coordinates corresponding to perturbations in co 2 , co 1 , 

and co 3 , respectively, from the equilibrium motion 


r 

7 

n 

CO 

4 > 


= Em/M 

= the nutation angle, i.e. the angle between the b 2 axis and H 
= nominal main body spin rate 

= angular velocities about the x, y, z axes respectively (i = 1,2.3) 

= position angle of the projection of dm R on the plane of the undeflected 
rotor disc 


<p = nutation damper displacement angle 
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Superscript 


( )' refers to component in hr' axes system 

( )’ indicates differentiation with respect to time 

Subscript 

b £ , i refers to particular main body axes (i = 1, 2, 3) 

0 refers to equilibrium value 

ss refers to steady state 


viii 



ACKNOWLEDGMENT 


The research reported here was conducted at the NASA Goddard 
Space Flight Center , Stabilization and Control Branch where the 
author was participating as a Faculty Fellow in the NASA-ASEE 
1971 Summer Faculty Fellowship Program. My research col- 
league was Dr. Joseph V. Fedor of the Stabilization and Control 
Branch whose continued interest is gratefully acknowledged. In 
addition, the author would like to express his appreciation to 
Dr. Arun K. Guha, Postdoctoral Research Associate in the Stabili- 
zation and Control Branch for the many illuminating discussions. 
Thanks are also extended to Mr. Peter Fuechsel of The Johns 
Hopkins University, Applied Physics Laboratory for calculating 
and supplying the various moments of inertia for the SAS-A system , 
and to Messrs. Leonard Moriarity and Jim Andary of the Stabiliza- 
tion and Control Branch for the computer programming. 


ix 



NUTATIONAL STABILITY OF A DUAL-SPIN SATELLITE UNDER 
THE INFLUENCE OF APPLIED REACTION TORQUES 

I. INTRODUCTION 

NASA Goddard Space Flight Center is currently directing the design and 
development of a series of Small Astronomy Satellites (SAS) whose purpose is 
to support orbitally based experiments in modern astronomy. The first of these 
satellites, SAS-A, is a dual-spin spacecraft which has the capability of scanning 
the entire celestial sphere to determine the relative position and intensity of 
X-ray emitting sources with respect to the fixed position of the stars. It is 
important that the attitude of the satellite be precisely known and maintained in 
order to accurately determine the location of the X-ray emitting sources. The 
attitude control system for the SAS-A satellite has been designed and developed 
by the Applied Physics Laboratory of the Johns Hopkins University. 

An analysis of the attitude motion and stability of such a dual-spin space- 
craft with damping only on the slowly spinning main part was reported previously . 1 
The resulting differential equations of rotational motion when linearized were an 
autonomous set of fifth order equations with constant coefficients. Analytical 
stability criteria were developed from these equations using the method of 
Routh-Hurwitz. 

Subsequent to this analysis, it was demonstrated by static and dynamic tests 
of the SAS-A momentum wheel that there was some energy dissipation in the 
shaft-momentum wheel assembly. A recent investigation 2 incorporated the 
effects of momentum wheel damping into the rotational equations of motion for 
the SAS-A spacecraft, and analytically considered the attitude stability of such a 
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system. Wheel energy dissipation was modeled by assuming the wheel could 
flex with two degrees of freedom relative to the hub. The nonlinear attitude 
equations derived for small wheel flexural motion were a ninth order nonautono- 
mous set. An averaging process was used to determine the zeroth and first 
order secular perturbations on the behavior of the system nutation angle. From 
this a general analytic stability criterion was established. A numerical evalua- 
tion of this criterion using parameters and measured wheel damping data for the 
Small Astronomy-A satellite indicated that stability about a zero degree nutation 
angle is insured by a factor of 128 under normal operating conditions. Numeri- 
cal integration of the nonlinear equations confirmed the analytic results for 
special cases. 

The SAS-A satellite was launched in December 1970. After initial despin 
maneuvers and a command to pulse the nutation damper magnet to a more 
optimum damping constant, the nutation angle amplitude was reduced from about 
1 degree to approximately 0.10-0.15 degree. Although this is within the mission 
requirement of 0.20 degree, it is about an order of magnitude greater than that 
predicted by prelaunch computer simulations. 3 Several hypotheses which might 
explain this discrepancy are currently under investigation. Among these are: 
possible periodic fluctuations in solar panel lengths due to thermal expansion/ 
contraction as the satellite rotates in the sunlight and which could excite nuta- 
tions; excessive mass imbalance of the satellite in-orbit configuration, in which 
the state of zero nutation angle may not correspond to the lowest energy state, 
and where nutation damping may approach a steady-state condition of minimum 
energy but with non-zero nutation; and thirdly, since the solar paddle attachments 
to the hub of the spacecraft are not rigidly locked, but rather spring-connected, 
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the effect of the (loose) solar paddles could be replaced by a constant reaction 
torque acting on the spacecraft hub. This could cause the satellite to reach an 
equilibrium motion about an axis slightly displaced from the nominal figure or 
symmetry axis , thus possibly explaining the anomalous nutational behavior of 
the SAS-A system. 

This third theory will be the basis for the investigation considered herein. 
Although the results will be applied directly to the SAS-A spacecraft, it should 
also be of interest to the designers of other dual-spin systems having non- 
rigidly connected solar paddles. 

II. ANALYSIS 

It is a well-known textbook problem that a rotating rectangular plate with a 
constant angular velocity can be forced to rotate about a nonprincipal axis of in- 
ertia under the influence of a reaction torque which acts about an axis perpen- 
dicular to the plate and rotates with it. 4 In 1963 Garber studied the stability of 
a gravity-gradient stabilized satellite under the influence of a constant disturbing 
(reaction) torque. 5 He concluded that if this torque were sufficient to produce a 
non-zero steady-state pitch (in-plane) motion, it could lead to instabilities in 
the roll-yaw out-of-plane motion. 5 

This idea may be applied to the motion of a spinning satellite with loosely 
attached solar paddles; a plausible model is shown in Figure 1. A possible equi- 
librium configuration is described by two body-axis components of the spacecraft 
spin vector that are non-zero: 
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<jl >2 - fi cos a 

0J t = o ( 1 ) 

a> 3 - 0. sin a 

where a is the angle between the axis of rotation and the desired system spin 
axis, and ft is the nominal value of the satellite spin rate. 

The equations of motion for an undamped dual-spin satellite system where 
the desired axis of rotation is the "y" axis can be represented by: 

Bd) 2 + cv 1 a) 3 (A - C) = L dz 

A Wj + o> 2 w 3 (C ~ B) - I R ^ s = L dj L (2) 

Cm 3 + <u) 2 (B “ A) + o) x I R 2 s = L dj 

If we calculate the disturbance torque, L d , necessary to give this system a 
forced rotation represented by Equation (1), we find that: 



u > 2 ~ ft cos a + S 

*>1 = € r ( 4 ) 

co 3 - ft sin a + i// 
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and substituting, Equation (4) back into the equations of motion (2). After 
linearization, the equations are of the form: 


where 


8 + Cj e = 0 

e + C 2 0 + C 3 5 = 0 ?■ 


ip + C 4 e - 0 


C x = (A “ C) n sin a/B 
c 2 = - [(B - C) n cos a. + I Rj s]/ A 
C 3 = - [(B - C) n sin a] /A 
C 4 = [(B - A) fi cos a + I R s]/C 


(5) 


If the first of Equations (5) is multiplied by C 4 , the third by -C x , and the 
result added, the linear dependence between 8 and i/' can be established. 


C 4 8 - C x ip = 0 


(6) 


or 


C 4 8 = C x ip + C 4 8 (0) ~C 1 <//(0) 


(?) 


where 8 (0), <//( 0) represent the initial conditions in the variational coordinates. 

After elimination of 8 as a dependent coordinate, the independent variational 

equations of motion can be written: 

. (C 2 C, * c 3 c^ c 3 'I 

e l F ) •]> = - F [C 4 S (0) - Cji/zCO)] 

^4 4 y 

I (8) 


i// + C 4 6 - 0 
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The homogeneous system has the following characteristic equation in the 
variable, 

- (c 2 c 4 + qc 3 ) = o 0) 

For a boundedness of the motion the roots of Equation (9) must be purely 
imaginary, or 

C 2 C 4 + C X C 3 < 0 (10) 

Inequality (10) may be expanded to take the form: 

B [(A ~ C) n sin a] [(§ - C) fi sin a] 

(U) 

+ C [(B - C) n cos a + I R2 s] [(B — A) Q cos a + I s] > 0 

Inequality (11) is obviously satisfied for the case of a SAS-type spacecraft where 
s »> Cl and approximately 90 percent of the system momentum is contributed 
by the rotor. 

For non-zero initial conditions in S or 0 it is seen from Equation (8) that 
non-zero steady-state values of S and <// would result, but also that e ss = o. It 
is clear that, even in the presence of damping, this system would never be 
asymptotically stable about the assumed equilibrium motion of Equation (1). Or, 
in other words , for small initial perturbations in the variational coordinates , 
regardless of the time elapsed, the system motion would never return exactly 
to its original equilibrium state. 

Up to this point the analysis has considered the forced rotation of a general 
asymmetric dual-spin satellite without energy dissipation on either the main 
part or the rotor. To consider the motion and stability of this system with 
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energy dissipation, first on the main part, and later on both parts of the space- 
craft, a specific model for the energy dissipation must be selected. For this 
purpose the nutation damping system employed in the SAS-type satellites will 
be considered. 


The satellite, wheel, and damping systems are illustrated in Figure 2. The 
elements of the attitude control system are: (1) the main part of the spacecraft, 
essentially a right circular cylinder where the nominal spin axis is the b 2 body 
axis, (2) a small momentum wheel or rotor assumed to be connected to the 
primary part near its center of mass, and whose spin axis is nominally parallel 
to the body axis, and (3) a pendulous-type nutation damper which is attached 
to the main part and constrained to move in a plane a distance t above the b x , b 3 
plane (determined by the body axes perpendicular to the symmetry axis, b 2 ). 

The damper is hinged about a torsion wire support which offers a restoring 
torque in addition to the dissipative torque associated with the damper motion. 


The equations of motion for this system under the reaction torque described 
above and with energy dissipation on the main body may be expressed: 


B ^ oj^co 3 (A C) L f 
Awj + co 2 u 3 (C - B) - w 3 : r 2 s 
Cu 3 + co 1 co 2 (B - A) + o) 1 I R2 s 



+ L 




( 12 ) 


in addition to an equation describing the single-degree-of-freedom damper. 
The damping terms, L f , and the damper equation as previously derived \ can 
be expanded to yield the following system of first order nonlinear differential 
equations under the assumptions that: 
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co i /co 2 « 1, 4> x « 1, 0J « 1, 

B w 2 + oi 1 w 3 (A - C) + m Tj (r j + r Q ) = 0 

A Wj + co 2 oj 3 (C “ B) - oj 3 I R2 s - 2m oj 2 r j T = L dj 
C oj 3 + WjOJj (B ~ A) + I Rs s - m rjT^j + mw ? 2 ^ = 

( ' m r j \ 

r° + —=—10! 

+ maJ 2 r 1 Fw 1 = - kc^j-K^j 


> ( 13 ) 


where the super-barred composite inertia terms now include the contributions 
of the main body, rotor, and the nutation damper. 

The linearized variational equations can again be obtained by substituting 
Equations (4) into (13) and linearizing, assuming: S « 1, e « l, 0 « i. The 

following set of linearized equations result: 


5 + Cj e + C = 0 

6 + C 2 t p + C 3 S + C 6 4> x = 0 

+ C 4 e + C + C 8 4> 1 = 0 


•n 


> ( 14 ) 


4> l + CgcJbj + C lo 0j + C n S + C 12 <// + C J3 e - 0 J 
where C x - C 4 have already been defined, after Equation (5), and, in addition: 


C 5 “ mr l( r l + r o ) /B 
C 6 = - 2m (fl cos a) r x r/A 

C 7 = - m r x r/C 
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Cg = mr 1 (fl 2 COS 2 a) r/C 
C 9 = k/mr 2 

C 1Q = K/ m r 2 + fi 2 cos 2 a (r Q /r 1 + m/M) 

C u = ( r o + r i)/ r i 

Cj 2 - - r/r, 

C 13 = r (Q cos a)/ r t 

Again after some algebraic manipulation to eliminate the dependent degree 
of freedom, it can be shown that the system characteristic equation corresponding 
to the homogeneous linear system can be expressed as: 

f (\) = a Q \ 4 + aj\ 3 + a 2 ^ 2 + a 3 \ + a 4 = 0 (lb) 

where 

a 0 _ 1 ^ 5^11 ” ^" 7^12 

a i = c 9 

a 2 = -(C 2 C 4 +c 1 c 3 )(i - C S C 1X -c 7 c 12 ) 

+ (^10 - ^8 ^ 12 ) 

- (-c 2 c 7 - c 3 c 5 + C 6 )(C 13 - c 1 c 11 - c 4 c 12 ) 

a 3 _ - Cg (C 2 C 4 + CjC 3 ) 

a 4 “ ^"2 ( _ ^"4 ^10 + ^ 8 ^ 13 ) 

A necessary condition for stability is that all of the coefficients of the 
characteristic equation be non-zero and have the same sign. After substituting 
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for the C's, it is found that 


a 0 = 1 " (m(r 0 + r^/B) - m r 2 /C 

and that for the actual SAS-type system the last two terms are much less 
than one. Also, 


a i - C 9 - k / m r 2 ; a j > 0 
if nutation damping is present. 

= C 9 (-C,C, -c,c,); 

if C 9 > 0, the conation that: -C 2 C„ -C, C, > 0 is the same as inequality 
(10) or (11) already discussed. 

For a « 1, but not zero, the necessary condition that a 4 > 0 can be re- 
duced, after some algebra, to: 

— +4+ K mr 2 m F(r 0 + r t ) (A - C)a 

r i M B-A+I^s/ft B [B - A + I s/fi] (16) 

When a — 0, = 0, and condition (16) is identical to the necessary and 

sufficient stability condition (18) of Reference 1. 

The condition that a 2 > 0 is automatically satisfied for the case when a = 0, 
and a 0 > o. For a « 1, it can be shown that a 2 > 0 when s »> ft. For a 
not restricted to be small the conditions that: a 2 > 0, a 3 > 0, a 4 > 0 are 
not as apparent, but could be developed and evaluated numerically. 

Wheel energy dissipation has been previously modeled by assuming that the 
momentum wheel could have two degrees of flexural freedom relative to the 
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main spacecraft hub. 2 Figure 3 describes the momentum wheel deflection 
geometry for the assumed small deflections a z and a x showing the displace- 
ment of a differential mass dm R on the deflected wheel. If it is assumed that the 
rotor is attached at or very near to the system center of mass and that the 
wheel is homogeneous (such that Ir 2 = 2I Rl ), together with the previous assump- 
tions on the magnitudes of the a 4> l and 4> x , the first order nonlinear equations 
of motion are obtained 2 : 

B u> 2 + co ^ co 3 ( A - C) + + r Q ) = L 2 

Awj + co 2 co 3 (C - B) - co 3 I R s - 2ma) 2 r 1 F ^ 

+ I R 1 {IX + («2 + s ) 2a z] sin st “ X + ("2 + s ) 2a x] cos st ! 

C 6J 3 + C0 1 CO 2 (B - A) + Wj I R2 s - + m w 2 2 r j T 

+ * R i { [® z + ("2 + s ) 2 a z] cos st + (X + ( w 2 + s ) 2a x] s ^ n st } 


(17) 

(18) 
} = Li 

(19) 

= l 3 


m r 


2 

1 


!r - mrj T co 3 + mr 1 (r 0 + r t ) co 2 + co 2 mr x 

M 



( 20 ) 


+ m aij fj T Wj 


- - K 4>\ 
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I R 1 a z + k R ®x + [ K r + I Rl 0 2 + s) 2 ] a z 

( 21 ) 

+ 1 { ["l ~ ("2 + 2s ) w 3 ] sin st + [ ^3 + (<U 2 + 2 s) WjJcos St) = o 

I R 1 -x + k R a, + [\ + I 8l 0 2 + S) 2 ] a x 

(22) 

+ I Rj { [ w 3 + ( w 2 + 2s ) ^i] sin st - [wj - (a> 2 + 2s) w 3 J cos st} =0 

If it is assumed that L 2 = L 3 = 0 and that L, = L di , Equations (17)-(22) may 

be linearized about the equilibrium motion: 

U) 2 - Q cos a + S, Wj = e, w 3 = Q sin a + \p 

$ l = a z ~ a x = 0 

to yield in terms of the previously defined constants: 

S + Cje + Cj.0! = 0 (23) 

e + C 2 <// + C 3 S + C 6 4> l + ^I Rl / [sin s t [ a z + (fi cos a + s) 2 a z ] 

(24) 

- cos st [a x + (fi cos a + s) 2 a x ]J = 0 

0 + C 4 e + C 7 + C g 4> x + ( I Ri /C){ cos s t [a, + (fi cos a + s) 2 a t ] 

J (25) 

+ sin s t [a x + (0 cos a + s) 2 a x ] } = 0 
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0 


(26) 


\ + C 9 ^ + Cj^j + C U S + C 12 0 + C 13 e 
+ ( k R /I Rl ) <*, + [K R /I Ri + (fi cos a + s) 2 ] a z 

+ [e - (Q cos a + 2 s )*/> + (fl sin a) S] sin st 
+ [\p + (fi cos a + 2s) e] cos st - (sin st)fi sin a (fl cos a + 2s) (27) 

+ + + cos a + s ) 2] a x 

+ [0 + (0 cos a + 2 s) e] sin st - [e - (fi cos a + 2s)\p + (fi sin a.) 5] cos st 

= - (cos st) fl sin a(f) cos a + 2s) (28) 

Equations (23)-(28) represent a ninth order non-homogeneous set of nonautono- 
mous differential equations. The stability of the homogeneous system could be 
analyzed using Floquet theory similar to the treatment of Mingori. 6 The appli- 
cation of Floquet theory for this problem would necessitate the use of a digital 
computer to study the stability point by point by varying different satellite and 
rotor inertia, spin, and damping system parameters in a systematic manner. 

Instead of the Floquet analysis the Lyapunov Reducibility Theorem may be 
applied to reduce the linear system of differential equations with periodic coeffi- 
cients to a kinematically similar system of autonomous differential equations 7 • 

8 

This method has been recently employed by Guha in studying the attitude sta- 
bility of dissipative dual-spin spacecraft. 
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The homogeneous system of linear nonautonomous differential equations 
can be expressed in standard state vector form: 

x = A( t ) x (29) 

where, for this system, A (t) represents a 12 x 12 matrix. According to the 
Lyapunov Reducibility Theorem 7,8 there is a bounded nonsingular transforma- 
tion matrix S (t) with a bounded inverse S ' *(t) such that a matrix, D, defined by 
the equation: 

D = S ' 1 AS - S ' (30) 

is a constant matrix. We can then relate, 


x = S(t)y (31) 

and obtain the system, 

y = Dy (32) 

The system of autonomous equations (32), is said to be kinematically similar to 
the original system, (29) and has similar stability properties since the charac- 
teristic roots of the matrix D are the characteristic exponents of the original 
equation, (29). According to a theorem of Lyapunov, every equation with periodic 

V. 

coefficients is reducible, and the major problem is to determine the bounded 
transformation, S(t). 

Guha has shown 8 in applying this technique to the system considered in 
Reference 2, that a suitable transformation of the form of Equation (31) can be 
expressed: 
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It can be verified that this same type of nonsingular bounded transformation can 
be applied to the present system, Equations (23)-(28). The following autonomous, 
and kinematically similar system in the variational coordinates S, e, 0, 4> v 


and /3 2 , results: 

S + Cj e + C s = 0 (34) 

e + C 2 0 + C 3 S + + C 14 ^ + C ls /3 2 + C 16 /3 2 = 0 (35) 

4 > +-C 4 e i +C 7 ^ 1 + Cg^ + C 17 /S 1 + C 18 /3j + C 19 /? 2 = 0 (36) 

A +C A +C 10^ + Cn s +C 12 0 + C 13 e = 0 (37) 

A + C 20 A + C 21 /3j + C 22/ S 2 + C 23 /3 2 + * + C 24 e = 0 (38) 


Z? 2 + C 20 /3 2 + C 21 /3 2 - C 22 ^ t - C 23 fi 1 + k - C 24 0 + C 25 s - n sin a (Cl cos a + 2s)(39) 
where C x -C have been defined previously and, in addition, 

Ci 4 - “2 I Ri s/5 

C 15 = 
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C 16 - I Ri Cl cos a (Cl cos a + 2 s) /A 

C 17 = I Bl /C 

C 18 = I Ri Cl cos a(fi cos a + 2s)/C 

C !9 = ^s/C 

^20 “ ^r/ I R 1 

C 2 i = (K R /I Ri ) + Q cos a (Cl + 2 s) 

C 22 = 2s 

^23 “ ^R S /^R! 

C 24 = Cl cos a + 2 s 

C 2S = Cl sin a 

Equation (38) is obtained by multiplying (27) by cos st and (28) by sin st and 
adding; Equation (39) is similarly obtained by multiplying (27) by sin st and (28) 
by -cos st and adding. 

The stability analysis of this system is a standard exercise in linear systems 
theory. The stability of the zero solution of the corresponding homogeneous 
system may be analyzed by Routh's method. The characteristic polynomial may 
be evaluated, and the system stability dependent on the roots of the character- 
istic equation all having negative real parts. 
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For the non-homogeneous system, the steady state solution may be obtained 
by assuming: 

S ss = K l- e ss = K 2- Si'ss = K 3- Ks = K 4’ = K S 

and B = K and substituting back into Equations (34)- (39). It can be seen that: 

r 2 ss 6 

^ss = C 24 C 2s/( C 25 + C 24 ’ 

Si 

'Pss ~ ' C 2 S ss 


and e =6 = B, = / 3 „ =0. An examination of the constants involved 

SS ^l s s ^1 S S r2 ss 

in S and \b indicates that the steady state motion of the kinematically similar 

S S ' s s 

system is bounded. 


III. NUMERICAL RESULTS 

The characteristic equation was developed for the linear autonomous 
homogeneous system corresponding to Equations (34) -(39). Equilibrium values 
of a 0 corresponding to the reportedly observed nutation angles of 0.1 degree 
and 0.2 degree were calculated as follows: 

C fi sin a n 

t a n y = t— 2 (4( 

tan y o I s + B Q cos a Q 

where the momentum components are approximated by: 

H 2 % Ir 2 s + B (° cos a ) 

Hj % 0 

H 3 ~ C (Q sin a) 
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since A ~ A, B ~ B, and C ~ C. 

Two cases of interest were considered, where, first one pair of SAS-A solar 
paddles was assumed to be loosely attached to the hub with the other pair rigidly 
connected, and secondly, vice versa. The moments of inertia for the two different 
cases are: 


Case I Case II 

A = 28.11929 kg-m 2 A = 15.18496 kg-m 2 

B = 17.747422 kg-m 2 B = 17.74742 kg-m 2 

C = 15.32054 kg-m 2 C = 28.634496 kg-m 2 

It is seen that from Equation (40), for Case I, a Q = 1°55' and for Case II, a 0 = 
l°r, for y 0 = 0.1 degree; also for y 0 = 0.2 degree, a 0 = 2°6\ and 3°49' for 
Casellandl respectively. 

The roots of the characteristic equation were calculated for the two cases 
using the following nominal SAS-A orbital parameters: 


T = 0.453 m. 
r Q = 0.0189 m. 
rj = 0.1735 m. 
m = 0.2451 kg. 

I Rz = 0.011525 kg-m 2 
I Rj = 0.0050799 kg-m 2 = I Rj 
s = 209.43 rad/sec. 
ft = 0.008727 rad/sec. 


k = 5.7 x 10“ 5 nt-m-sec/rad. 

K = 5.8 x 10" 5 nt-m/rad. 

*k R = 6.779 x 10 -3 nt-m-sec/rad. 
*K r = 71.61074 nt-m/rad. 

M = 114 kg. 


*Determined experimentally at The Johns Hopkins University, Applied Physics Laboratory (see 
Ref. 2). 
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For this system, after elimination of the dependence between coordinates, 
similar to Equations (6) and (7), an eighth order characteristic equation results. 
The roots were obtained using digital computations with a polynomial root-finder 
subroutine on the XDS-9300 computer. The results are summarized as follows 
for y Q = 0.1 degree (for y 0 = 0.2 degree, changes are observed only in the 
third or fourth significant figure): 


Real Part 

Imaginary Part - (Rad/sec) 

Case I 

-.37874 x 10' 2 

±.88380 x 10' 1 

-.8058 x 10' 4 

±.11455 

-.58874 

±.27962 x 10 2 

-.58843 

±.44670 x 10 3 

Case II 

-.38139 x 10' 2 

±.88468 x 10' 1 

-.4834 x 10' 4 

±.11374 

-.58875 

±.27963 x 10 2 

-.58844 

±.4467 x 10 3 


It is seen that all of the roots for both cases contain only negative real parts. 
The higher frequency modes are essentially the same for the two cases. The 
time constant for the least damped mode is about 12,400 seconds (3.44 hrs.) for 
the first case and 20,700 seconds (5.75 hrs.) for the second case. This repre- 
sents a considerable degradation in the performance of the nutation damping 
system about the equilibrium position y 0 = 0 previously considered. 1 • 2 
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The highest frequency mode for the kinematically similar system has the 
frequency of 446.7 rad/sec, approximately 2.14 times the wheel relative angular 
velocity of 209.4 rad/sec. It should be noted that in the original system the 
wheel deflection equations, Equations (27) and (28), have forcing terms on the 
right hand side at frequency s = 209.4 rad/sec. The possibility of forced 
resonance in the original system might exist, especially as the highest frequency 
of the kinematically similar system approaches 2s. This eventuality could be 
examined more closely by numerical integration of the system equations. 

It has been noted in the SAS-A post launch performance, that in addition to 
the somewhat larger nutation angles experienced than predicted by pre-launch 
analyses, a noticeable degradation of the nutation damping time constant occurs. 
Time constants as large as three hours or more have been observed. 9 While 
the present analysis does not claim an exact correlation with the anomalous 
observed behavior, it does suggest a possible mechanism for future problems, 
especially in connection with high pointing accuracy dual-spin missions. 

It has also been demonstrated by computer simulations that external torques 
having amplitudes of about 4000 dyne-cm would be necessary to produce nuta- 
tional motion of the SAS-A satellite with amplitudes of 0.1 degree. 9 Torques of 
this magnitude would be difficult to justify using even conservative aerodynamic 
density values at the SAS-A orbital altitude. The magnitude of the reaction 
torque, L dj , as given by Equation (3) was calculated for the two cases previously 
considered; for Case I, = -7090 dyne-cm, whereas for Case II, L di = -3620 
dyne-cm. There is a strong suggestion, therefore, that the solar panels and not 
the atmosphere are providing the higher amplitude torques associated with the 
anomalous nutational motion. 
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IV. CONCLUSIONS — CLOSING COMMENTS 


Based on the present stability analysis and numerical results the following 
concluding comments can be made: 

1. The motion of a dual-spin satellite system, under the influence of 
applied reaction torques , about an axis of rotation slightly displaced 
from the nominal spin axis remains bounded. For small initial per- 
turbations, the system, however, is not asymptotically stable about the 
original equilibrium state. 

2. For the case of energy dissipation only on the main body it is seen that 
the necessary Routh stability criteria for this motion are satisfied; for 
the general case of damping on both the main spacecraft and the wheel, 
stability is indicated by a numerical examination of the roots of the 
system characteristic equation of the kinematically similar system. 

3. The Lyapunov Reducibility Theorem can be used in studying the stability 
of the original nonautonomous system for the case of energy dissipation 
on both parts of the dual-spin system. 

4. A possible mechanism for the somewhat anomalous nutational motion 
and large nutational time constants observed in the post-launch SAS-A 
flight data is suggested. For this system it is assumed that the applied 
reaction torques would be provided by either pair of non-rigidly attached 
solar paddles. 
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Figure 2. Elements of SAS-A Attitude Control System 
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gure 3. Momentum Wheel Deflection Geometry 


